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Abstract. We discuss the production of charmonium in nuclear collisions within the 
i framework of the statistical hadronization model. We demonstrate that the model 

C*~) | reproduces very well the availble data at RHIC. We provide predictions for the LHC 

energy where, dependently on the charm production cross section, a dramatically 
different behaviour of charmonium production as a function of centrality might be 
expected. We extend our predictions for charm production towards the threshold 
Q ■ energies, where charm is expected to be measured at the future FAIR facility. 

a. 

Charmonium production is considered, since the original proposal more than 20 
years ago about its suppression in a Quark-Gluon Plasma (QGP) pQ, an important 
probe to determine the degree of deconfinement reached in the fireball produced in 
ultra-relativistic nucleus-nucleus collisions. In the original scenario of J /if) suppression 
via Debye screening [I] it is assumed that the charmonia are rapidly formed in initial 
hard collisions but are subsequently destroyed in the QGP (see an update in ref. [2]). 

In recent publications [3] we have demonstrated that the data on J/if) and if)' 
production in nucleus- nucleus collisions at the SPS (y/SNN ~ 17 GeV) and RHIC 
(v^viv = 200 GeV) energies can be well described within the statistical hadronization 
^ ■ model proposed in [4j and have provided predictions for the LHC energy (a/saT/v=5.5 
TeV) and for energies close to threshold {^s^n ~ 6 GeV). 

In our statistical hadronization model (SHM) [U [5l [3] we assume that the charm 
quarks are produced in primary hard collisions and that their total number stays 
constant until hadronization at chemical freeze-out. Another important element 
is thermal equilibration in the QGP, at least near the critical temperature, T c . 
Measurements of elliptic flow of J /if) will be a crucial check of this scenario. 

The model has the following input parameters: i) the charm production cross 
section in pp collisions, taken either from NLO pQCD calculations [HI E] or from 
experiment [8]; ii) characteristics at chemical freeze-out: temperature, T, baryochemical 
potential, /if,, and volume corresponding to one unit of rapidity Va v =i, extracted from 
thermal fits of non-charmed hadrons [9]. 

In Fig. CD we present the rapidity dependence of the nuclear modification factor R J Jj/- 
While earlier [3] we have compared data [10] to our model predictions for the pQCD 
charm production cross section (6] (dashed line in Fig. [1]), we use here as alternative 
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the charm cross section as measured by PHENIX in pp collisions [8] and consider in 
addition shadowing for Au-Au collsions as extracted from recent dAu data [11] (assuming 
teh deviation of R J Ja u from unity is due to shadowing). Also in this case, our model 
describes the observed suppression and its rapidity dependence. The maximum of R'aa 
at midrapidity is in our model due to the enhanced generation of charmonium around 
mid-rapidity, determined by the rapidity dependence of the charm production cross 
section. In this sense, the above result constitutes strong evidence for the statistical 
generation of J/ip at chemical freeze-out. 
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Figure 1. Rapidity dependence of R J Jt for two centrality classes. The data from the 
PHENIX experiment [TU] (symbols with errors) are compared to calculations (lines, 
see text) . The shaded area corresponds to calculations for the lower limit of the charm 
cross section as measured by PHENIX [SJ, with our shadowing scenario. 

The centrality dependence of R'aa a ^ midrapidity is shown in the left panel of 
Fig. [2J Our calculations approach the value in pp collisions around N part =50, which 
corresponds to an assumed minimal volume for the creation of QGP of 400 fm 3 [3]. 
The model reproduces very well the decreasing trend versus centrality seen in the RHIC 
data [TD] . Note that, in our model, the centrality dependence of the nuclear modification 
factor arises entirely as a consequence of the still rather moderate rapidity density of 
initially produced charm quark pairs at RHIC (dN c5 /dy=1.6 for central collisions, using 
the FONLL charm production cross section [6]). At the much higher LHC energy the 
charm production cross section (including shadowing in PbPb collisions [7]) is expected 
to be about an order of magnitude larger. As a result, the opposite trend as a function of 
centrality is predicted, with R^a exceeding unity for central collisions. A significantly 
larger enhancement of about a factor of 2 is obtained if the charm production cross 
section is two times larger than presently assumed, as seen in the right panel of Fig. [2j 
where we show the J/ip production relative to the number of initially produced cc pairs. 

A comprehensive set of model predictions from RHIC energy down to the charm 
production threshold is presented in Fig.[3]for central Au-Au collisions (A^ art =350). The 
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Figure 2. Centrality dependence of R/a f° r R-HIC and LHC energies (left panel) 
and of the J/t/j rapidity density at LHC relative to the number of initially produced 
cc pairs (right panel, curves labelled by the da c5 /dy) at midrapidity. 



left panel shows our predictions for the energy dependence of midrapidity yields (relative 
to the cc yield) for various charmed hadrons. The most striking behavior is observed 
for the production of A+ baryons: their yield rises significantly towards lower energies. 
In our approach this is caused by the increase in baryochemical potential towards lower 
energies (coupled with the charm neutrality condition). A similar behavior is seen for the 
baryon. The relative production yields of D-mesons depend on their quark content 
and depend on energy only around threshold. These results emphasize the importance 
of measuring, at low energies, in addition to D-mesons, also the yield of charmed baryons 
to get a complete measure of the total charm production cross section. 

One of the motivations for the study of charm production at low energies 
was the expectation [121 13] to provide, by a measurement of D-meson production 
near threshold, information on their possible in-medium modification near the phase 
boundary. However, the cross section a c5 is governed by the mass of the charm quark 
m c pa 1.3 GeV, which is much larger than any soft Quantum Chromodynamics (QCD) 
scale such as Aqcd- Therefore we expect no medium effects on this quantity. The 
much later formed D-mesons, or other charmed hadrons, may well change their mass in 
the hot medium. Whatever the medium effects may be, they can, because of the charm 
conservation, a c5 = ^(o~d + o~a c + o~s c + ...) + (a, rjc + aj/ 1 p + a Xc + ...), in first order only lead 
to a redistribution of charm quarks [3] . This argument is essentially model- independent 
and applies equally at all energies. This is demonstrated in the right panel of Fig. [31 
where we plot the relative change of the yields for different in-medium mass scenarios 
(see ref. [3] for details and references therein) compared to the case of vacuum masses. 
In contrast, the yields of charmonia do vary and this is more prominent at threshold 
energies (see Fig. [3]). 

We have shown that the statistical hadronization model describes well the measured 
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Figure 3. Left panel: energy dependence of the production yield relative to the 
number of cc pairs for open and hidden charm hadrons (for mesons the open symbols 
are for the antiparticles) . Right panel: relative change in the production yield of all 
open charm hadrons and of J/tp meson considering different scenarios for in-medium 
mass modifications (see text). 



decrease with centrality and the rapidity dependence of R A X at RHIC energy. 
Extrapolation to LHC energy leads, contrary to the observations at RHIC, to R~aa 
increasing with collision centrality. The increasing importance at lower energies of A c 
production was discussed and provides a challenge for future experiments. We have 
also shown that possible modifications of charmed hadrons in the hot hadronic medium 
do not lead to measurable changes in the cross sections for production of hadrons with 
open charm. A possible influence of medium effects can be seen, however, in the yields 
of charmonium. 
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